Evolutionary dynamics of hyperbolic language.

Models of evolution of simple languages have typically assumed full alignment of the speaker and listeners interests, with perfect understanding representing the optimal outcome for both parties. In more realistic settings, communicating individuals will often desire different outcomes from one another. Previous work has shown that misalignment of speaker-listener interests reduces the maximum informativeness among Nash-equilibrium languages, and that multiple equilibrium languages (with different degrees of informativeness) are supported. We study the stochastic evolutionary dynamics of signaling games in which the alignment of speaker-listener interests can vary. We find that increased misalignment of speaker-listener interests is associated with a decrease in information transmission. Moreover, the most common languages to evolve are typically the most informative languages supportable as static Nash equilibria, suggesting a solution to the 'equilibrium selection problem'. In addition, our dynamics reveal the mechanism by which less informative languages evolve: words that previously signaled intense states come to be used hyperbolically for less intense states, with listeners' interpretation of these newly-ambiguous words evolving downward in response. We ground our results in linguistic data on intensifiers such as so and very, words which have unique dynamics-with constant recycling and innovation that match our theoretical results well.

Longitudinal study of 2 patients with cyclic thrombocytopenia, STAT3 and MPL mutations.

Cyclic thrombocytopenia (CTP) is a rare disease of periodic platelet count oscillations. The pathogenesis of CTP remains elusive. To study the underlying pathophysiology and genetic and cellular associations with CTP, we applied systems biology approaches to 2 patients with stable platelet cycling and reciprocal thrombopoietin (TPO) cycling at multiple time points through 2 cycles. Blood transcriptome analysis revealed cycling of platelet-specific genes, which are in parallel with and precede platelet count oscillation, indicating that cyclical platelet production leads platelet count cycling in both patients. Additionally, neutrophil and erythrocyte-specific genes also showed fluctuations correlating with platelet count changes, consistent with TPO effects on hematopoietic progenitors. Moreover, we found novel genetic associations with CTP. One patient had a novel germline heterozygous loss-of-function (LOF) thrombopoietin receptor (MPL) c.1210G>A mutation, and both had pathogenic somatic gain-of-function (GOF) variants in signal transducer and activator of transcription 3 (STAT3). In addition, both patients had clonal T-cell populations that remained stable throughout platelet count cycles. These mutations and clonal T cells may potentially involve in the pathogenic baseline in these patients, rendering exaggerated persistent thrombopoiesis oscillations of their intrinsic rhythm upon homeostatic perturbations. This work provides new insights into the pathophysiology of CTP and possible therapies.

A Blueprint for Identifying Phenotypes and Drug Targets in Complex Disorders with Empirical Dynamics.

A central challenge in medicine is translating from observational understanding to mechanistic understanding, where some observations are recognized as causes for the others. This can lead not only to new treatments and understanding, but also to recognition of novel phenotypes. Here, we apply a collection of mathematical techniques (empirical dynamics), which infer mechanistic networks in a model-free manner from longitudinal data, to hematopoiesis. Our study consists of three subjects with markers for cyclic thrombocytopenia, in which multiple cells and proteins undergo abnormal oscillations. One subject has atypical markers and may represent a rare phenotype. Our analyses support this contention, and also lend new evidence to a theory for the cause of this disorder. Simulations of an intervention yield encouraging results, even when applied to patient data outside our three subjects. These successes suggest that this blueprint has broader applicability in understanding and treating complex disorders.

Population structure across scales facilitates coexistence and spatial heterogeneity of antibiotic-resistant infections.

Antibiotic-resistant infections are a growing threat to human health, but basic features of the eco-evolutionary dynamics remain unexplained. Most prominently, there is no clear mechanism for the long-term coexistence of both drug-sensitive and resistant strains at intermediate levels, a ubiquitous pattern seen in surveillance data. Here we show that accounting for structured or spatially-heterogeneous host populations and variability in antibiotic consumption can lead to persistent coexistence over a wide range of treatment coverages, drug efficacies, costs of resistance, and mixing patterns. Moreover, this mechanism can explain other puzzling spatiotemporal features of drug-resistance epidemiology that have received less attention, such as large differences in the prevalence of resistance between geographical regions with similar antibiotic consumption or that neighbor one another. We find that the same amount of antibiotic use can lead to very different levels of resistance depending on how treatment is distributed in a transmission network. We also identify parameter regimes in which population structure alone cannot support coexistence, suggesting the need for other mechanisms to explain the epidemiology of antibiotic resistance. Our analysis identifies key features of host population structure that can be used to assess resistance risk and highlights the need to include spatial or demographic heterogeneity in models to guide resistance management.

Turbulent coherent structures and early life below the Kolmogorov scale.

Major evolutionary transitions, including the emergence of life, likely occurred in aqueous environments. While the role of water's chemistry in early life is well studied, the effects of water's ability to manipulate population structure are less clear. Population structure is known to be critical, as effective replicators must be insulated from parasites. Here, we propose that turbulent coherent structures, long-lasting flow patterns which trap particles, may serve many of the properties associated with compartments - collocalization, division, and merging - which are commonly thought to play a key role in the origins of life and other evolutionary transitions. We substantiate this idea by simulating multiple proposed metabolisms for early life in a simple model of a turbulent flow, and find that balancing the turnover times of biological particles and coherent structures can indeed enhance the likelihood of these metabolisms overcoming extinction either via parasitism or via a lack of metabolic support. Our results suggest that group selection models may be applicable with fewer physical and chemical constraints than previously thought, and apply much more widely in aqueous environments.

Tunable wrinkling of thin nematic liquid crystal elastomer sheets.

Instabilities in thin elastic sheets, such as wrinkles, are of broad interest both from a fundamental viewpoint and also because of their potential for engineering applications. Nematic liquid crystal elastomers offer a new form of control of these instabilities through direct coupling between microscopic degrees of freedom, resulting from orientational ordering of rodlike molecules, and macroscopic strain. By a standard method of dimensional reduction, we construct a plate theory for thin sheets of nematic elastomer. We then apply this theory to the study of the formation of wrinkles due to compression of a thin sheet of nematic liquid crystal elastomer atop an elastic or fluid substrate. We find the scaling of the wrinkle wavelength in terms of material parameters and the applied compression. The wavelength of the wrinkles is found to be nonmonotonic in the compressive strain due to the presence of the nematic. Finally, due to soft modes, the critical stress for the appearance of wrinkles can be much higher than in an isotropic elastomer and depends nontrivially on the manner in which the elastomer was prepared.

Large-scale quantitative profiling of the Old English verse tradition.

The corpus of Old English verse is an indispensable source for scholars of the Indo-European tradition, early Germanic culture and English literary history. Although it has been the focus of sustained literary scholarship for over two centuries, Old English poetry has not been subjected to corpus-wide computational profiling, in part because of the sparseness and extreme fragmentation of the surviving material. Here we report a detailed quantitative analysis of the whole corpus that considers a broad range of features reflective of sound, metre and diction. This integrated examination of fine-grained features enabled us to identify salient stylistic patterns, despite the inherent limitations of the corpus. In particular, we provide quantitative evidence consistent with the unitary authorship of Beowulf and the Cynewulfian authorship of Andreas, shedding light on two longstanding questions in Old English philology. Our results demonstrate the usefulness of high-dimensional stylometric profiling for fragmentary literary traditions and lay the foundation for future studies of the cultural evolution of English literature.

Effects of motion in structured populations.

In evolutionary processes, population structure has a substantial effect on natural selection. Here, we analyse how motion of individuals affects constant selection in structured populations. Motion is relevant because it leads to changes in the distribution of types as mutations march towards fixation or extinction. We describe motion as the swapping of individuals on graphs, and more generally as the shuffling of individuals between reproductive updates. Beginning with a one-dimensional graph, the cycle, we prove that motion suppresses natural selection for death-birth (DB) updating or for any process that combines birth-death (BD) and DB updating. If the rule is purely BD updating, no change in fixation probability appears in the presence of motion. We further investigate how motion affects evolution on the square lattice and weighted graphs. In the case of weighted graphs, we find that motion can be either an amplifier or a suppressor of natural selection. In some cases, whether it is one or the other can be a function of the relative reproductive rate, indicating that motion is a subtle and complex attribute of evolving populations. As a first step towards understanding less restricted types of motion in evolutionary graph theory, we consider a similar rule on dynamic graphs induced by a spatial flow and find qualitatively similar results, indicating that continuous motion also suppresses natural selection.

Interfacial fluid instabilities and Kapitsa pendula.

The onset and development of instabilities is one of the central problems in fluid mechanics. Here we develop a connection between instabilities of free fluid interfaces and inverted pendula. When acted upon solely by the gravitational force, the inverted pendulum is unstable. This position can be stabilized by the Kapitsa phenomenon, in which high-frequency low-amplitude vertical vibrations of the base creates a fictitious force which opposes the gravitational force. By transforming the dynamical equations governing a fluid interface into an appropriate pendulum-type equation, we demonstrate how stability can be induced in fluid systems by properly tuned vibrations. We construct a "dictionary"-type relationship between various pendula and the classical Rayleigh-Taylor, Kelvin-Helmholtz, Rayleigh-Plateau and the self-gravitational instabilities. This makes several results in control theory and dynamical systems directly applicable to the study of tunable fluid instabilities, where the critical wavelength depends on the external forces or the instability is suppressed entirely. We suggest some applications and instances of the effect ranging in scale from microns to the radius of a galaxy.

Microorganism billiards in closed plane curves.

Recent experiments have shown that many species of microorganisms leave a solid surface at a fixed angle determined by steric interactions and near-field hydrodynamics. This angle is completely independent of the incoming angle. For several collisions in a closed body this determines a unique type of billiard system, an aspecular billiard in which the outgoing angle is fixed for all collisions. We analyze such a system using numerical simulation of this billiard for varying tables and outgoing angles, and also utilize the theory of one-dimensional maps and wavefront dynamics. When applicable we cite results from and compare our system to similar billiard systems in the literature. We focus on examples from three broad classes: the ellipse, the Bunimovich billiards, and the Sinai billiards. The effect of a noisy outgoing angle is also discussed.

Microscale locomotion in a nematic liquid crystal.

Microorganisms often encounter anisotropy, for example in mucus and biofilms. We study how anisotropy and elasticity of the ambient fluid affects the speed of a swimming microorganism with a prescribed stroke. Motivated by recent experiments on swimming bacteria in anisotropic environments, we extend a classical model for swimming microorganisms, the Taylor swimming sheet, actuated by small-amplitude traveling waves in a three-dimensional nematic liquid crystal without twist. We calculate the swimming speed and entrained volumetric flux as a function of the swimmer's stroke properties as well as the elastic and rheological properties of the liquid crystal. These results are then compared to previous results on an analogous swimmer in a hexatic liquid crystal, indicating large differences in the cases of small Ericksen number and in a nematic fluid when the tumbling parameter is near the transition to a shear-aligning nematic. We also propose a novel method of swimming or pumping in a nematic fluid by passing a traveling wave of director oscillation along a rigid wall.

Minimal model for transient swimming in a liquid crystal.

When a microorganism begins swimming from rest in a Newtonian fluid such as water, it rapidly attains its steady-state swimming speed since changes in the velocity field spread quickly when the Reynolds number is small. However, swimming microorganisms are commonly found or studied in complex fluids. Because these fluids have long relaxation times, the time to attain the steady-state swimming speed can also be long. In this article we study the swimming startup problem in the simplest liquid crystalline fluid: a two-dimensional hexatic liquid crystal film. We study the dependence of startup time on anchoring strength and Ericksen number, which is the ratio of viscous to elastic stresses. For strong anchoring, the fluid flow starts up immediately but the liquid crystal field and swimming velocity attain their sinusoidal steady-state values after a time proportional to the relaxation time of the liquid crystal. When the Ericksen number is high, the behavior is the same as in the strong-anchoring case for any anchoring strength. We also find that the startup time increases with the ratio of the rotational viscosity to the shear viscosity, and then ultimately saturates once the rotational viscosity is much greater than the shear viscosity.

Locomotion and transport in a hexatic liquid crystal.

The swimming behavior of bacteria and other microorganisms is sensitive to the physical properties of the fluid in which they swim. Mucus, biofilms, and artificial liquid-crystalline solutions are all examples of fluids with some degree of anisotropy that are also commonly encountered by bacteria. In this article, we study how liquid-crystalline order affects the swimming behavior of a model swimmer. The swimmer is a one-dimensional version of G. I. Taylor's swimming sheet: an infinite line undulating with small-amplitude transverse or longitudinal traveling waves. The fluid is a two-dimensional hexatic liquid-crystalline film. We calculate the power dissipated, swimming speed, and flux of fluid entrained as a function of the swimmer's wave form as well as properties of the hexatic film, such as the rotational and shear viscosity, the Frank elastic constant, and the anchoring strength. The departure from isotropic behavior is greatest for large rotational viscosity and weak anchoring boundary conditions on the orientational order at the swimmer surface. We even find that if the rotational viscosity is large enough, the transverse-wave swimmer moves in the opposite direction relative to a swimmer in an isotropic fluid.